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The quantum criticality of the two-lead two-channel pseudogap Anderson model is studied. Based 
on the non-crossing approximation, we calculate both the linear and nonlinear conductance of the 
model at finite temperatures with a voltage bias and a power-law vanishing conduction electron 
density of states, Pc(w) oc |a; — (0 < r < 1) near the Fermi energy fip. Equilibrium and non¬ 

equilibrium quantum critical properties at the two-channel Kondo (2CK) to local moment (LM) 
phase transition are addressed by extracting universal scaling functions in both linear and non-linear 
conductances, respectively. Clear distinctions are found on the critical exponents between linear and 
non-linear conductance. The implications of these two distinct quantum critical properties for the 
non-equilibrium quantum criticality in general are discussed. 

PACS numbers: 


I. INTRODUCTION. 


Quantum phase transitions (QPTs)[lj, the zero- 
temperature phase transitions due to quantum fluctu¬ 
ations, are of fundamental importance in condensed 
matter systems. Near the transitions, these systems 
show non-Fermi liquid behaviors manifested in univer¬ 
sal power-law scaling in all thermodynamic observables. 
Recently, QPTs in quantum impurity problems^], have 
attracted much attention recently due to their relevance 
for the nano-systems, such as: quantum dots[3|, realized 
experimentally. The well-known Kondo effect]^, the an¬ 
tiferromagnetic spin correlations between impurity and 
conduction electrons, plays a crucial role in understand¬ 
ing their low temperature behaviors. New scaling laws 
are expected to occur when the QPTs are associated with 
the Kondo breakdown in these systems either in equilib¬ 
rium (at finite temperatures) or under non-equilibrium 
conditions (at finite voltage bia^. Of particular interest 
lies in QPTs out of equilibrium!^ where distinct univer¬ 
sal scalings are expected in contrast to the counterparts 
in equilibrium. 

A fascinating playground to address this issue is the ex¬ 
otic two-channel Kondo (2CK)[mil systems with non- 
Fermi liquid ground state due to overscreening of s = 1/2 
impurity spin by two independent conduction reservoirs. 
Much of the theoretical effort has been made for the 


2CK physics, including: via Bethe ansatz [12|, confor¬ 
mal field theory 0 , bosonization 0 and the numeri¬ 
cal renormalization group 0 . Experimentally, the 2CK 
ground state has been realized in semiconductor quan¬ 
tum dots Il6ll. magnetically doped nanowires, and metal¬ 
lic glasses [171 0. Recently, Kondo physics in magnet¬ 
ically doped graphene has attracted much attention for 
the possible 2CK physics as well as the pseudogap local 
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FIG. 1: Schematic 3D phase diagram of 2CK-LM quantum 
phase transition of the two-lead 2-channel pseudogap Ander¬ 
son model as functions of voltage bias V, temperature T and 
power-law exponent r of the pseudogap fermion bath. 


density of states (LDOS) Pc(w) which vanishes linearly 
due to the Dirac spectrum: pdto) oc |a;|'' with r = 1 |19- 
. This leads to QPT from the Kondo screened phase to 
the unscreened local moment (LM) phase with decreas¬ 
ing Kondo correlation due to insufficient DOS of con¬ 
duction electrons. In fact, such transitions exist in the 
more general framework of pseudogap Kondo (or Ander¬ 
son) models with power-law exponent 0 < r < 1, which 
have been extensively studied However, rela¬ 

tively less is known on the more exotic 2CK-LM QPT in 
the pseudogap Anderson (Kondo) models 27], [2^ , in par¬ 
ticular, when the system is subject to a non-equilibrium 
condition!^, 0. 
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0 < r < 1, where D = 1 is the bandwidth cutoff. 
Graphene and highTc-cuprate superconductors are pos¬ 
sible realizations of the pseudogap leads with r = 1, 
while semiconductors with soft gaps are candidates with 
0 < r < 1. The Hamiltonian reads [1, [s^ 

T,a,(7,k (T,r 

cr,(7' ,T,r' 

+ E {VkZctlda + h.C.) 

T,a,(7,k 

( 1 ) 


FIG. 2: The equilibrium impurity spectral function p{u}) (in 
arbitrary unit) versus energy u (in units of the bandwidth D) 
with 1/ = 0, r = 0.28, £d = -0.2, T = 5 X 10"^ (a). The 
spectral function p(a;) in the full energy range. For r < 

(re ~ 0.115), the Kondo peaks are located at w = 0 (Fermi 
energy), while for r > r^, dips near cu = 0 are developed, (b). 
Same plots as in (a) with energies close to the peaks/dips. 


Motivated by these developments, we address the 2CK- 
LM quantum phase transition in the two-channel pseu¬ 
dogap Anderson model both in and out of equilibrium 
in a Kondo quantum dot subject to a voltage bias and 
finite temperature. By studying experimentally accessi¬ 
ble steady-state transport, we search for universal scal¬ 
ings in both linear and non-linear conductance near crit¬ 
icality via a large-N method based on Non-Crossing 
Approximation (NCA) [^, 31-33|, a reliable approach 


for multi-channel Kondo systems with non-Fermi liq¬ 
uid ground states. A fundamentally important but less 
addressed issue-non-equilibrium quantum criticality- is 
emphasized here by identifying and comparing differ¬ 
ent universal scaling behaviors between equilibrium (zero 
bias) and non-equilibrium (finite bias) conductances near 
the 2CK-LM transition[33. 


II. THE MODEL HAMILTONIAN. 

The Hamiltonian of the two-lead two-channel single¬ 
impurity pseudogap Anderson model is formulated 
within the Non-Crossing Approximation (NCA) [i, lul¬ 
ls^, a large-N approach based on the SU{N) x SU{M) 
generalization of the SU (2) model with A —>■ oo, M —oo 
being the number of degenerate flavors of spins a =t 
N and the number of Kondo screening channels 
K = 1,2, •••M. In the physical SU[2) two-channel 
Kondo system, N = M = 2. The two leads are de¬ 
scribed by a power-law vanishing density of states (DOS) 
at Fermi energy defined as Pc(w) ^ |a;|’’0(ll — |a;|) with 


where iia=L/R = ^Vl2 is the chemical potential of lead 
a = L/R. The operators create (destroy) 

an electron in the leads with momentum k. Spin flavors 
are represented by cr, ct' = 1, • • • A and r, r' = 1, • • • M 
corresponds to M independent electron reservoirs. Here, 
A = M = 2. The U term describes the on-site Coulomb 
energy on the impurity, and represents the hybridiza¬ 
tion strength between the graphene electrons and the im¬ 
purity. 

The Hamiltonian Eq. m can be solved via NCA in the 
large-U limit 17 —> oo with the following Hamiltonian in 
the pseudofermion slave-boson representation: 

H = E “ fJ-aKlAr + Cd E 

kara. <7 

karex. 

( 2 ) 


where the local (impurity) electron operator rff ^ is de¬ 
composed in the pseudofermion representation as a prod¬ 
uct of pseudofermion /J and a slave-boson br- d\ ^ = 
f^br subject to the local constraint Q = fafer + 
ErKbr = 1 via the Lagrange multiplier A to ensure 
single occupancy on impurity. 

Here, we employ the NCA to address the equilib¬ 
rium and non-equilibrium transport at criticality based 
on Eq. O- This approach has been known to cor¬ 
rectly capture the non-Fermi liquid properties of the 
two-channel Anderson model[3, 28, 2^. Recently, it has 
been generalized to address the 2CK-LM crossover in 
non-equilibrium transport in a voltage-biased 2-channel 
pseudogap Anderson model with r = 1, relevant for 
graphene [30|. We generalize this approach here further to 
the voltage-biased 2-channel pseudogap Anderson model 
with 0 < r < 1. 

Within NCA, the Green’s functions for the conduction 
electrons pseudo-fermions Gfa{t), and the slave 


bosons D{t) are given bv|8l. l32l|: 


< {c^{t),cl{0)} >ix 



















FIG. 3: The non-equilibrium impurity spectral function 
(in arbitrary unit) with F = 0.3D, td = —0.3D, r = 0.05, 
V = 0.038Z), T = 5 X 10~^D. The Kondo peak splits into 
two peaks where the width of the two splited peaks equals to 
bias V (see Inset). 



FIG. 4: The equilibrium (linear) conductance G(0, T) (in unit 
of 2elh) versus temperature T (in unit of D = 1) for various 
r with F = 0.3D, td = —0.3D. 


( 3 ) 

G% = -* < /.(f)/](0) >.x,G% = i < fliO)Ut) >^x, 

( 4 ) 

and 

D> = -i< >.x,D< = i < bU0)br{t) >ix, 

( 5 ) 

where the notation < and > represents lesser and greater 
Green function. The lesser (greater) Green functions can 
be written as [1, : 

( 6 ) 


Here, rQ,(a;) = TaPc,a{uj) with r„ = 7r|I4“p with 
Pc.a(w) = -^ImGciuj) = ^j^\uj\^9{D -\uj- p,a\) and 
D = 1 being the bandwidth of the conduction bath, and 
/(<^) = ^^le Fermi function. 

The relation between greater Green function and re¬ 
tarded Green function are, D^{w) = 2iImD^{w), G^^ = 
2iImG'j^, where the self-energies, n^(r(;) = 2iImIV{w) 
and = 2z/mE^^(w). The NGA expressions for 

the self energies of retarded Greens functions for pseudo¬ 
fermion Gj^(w) = [w —Cd —E’'(w)]“^ and for slave-boson 
D'~{u;) = [tu — n’’(a;)]“^ are given by[l, 3^, 3^ : 

E’'(w) = -^ f dera{u;-e-Ha) f{e-uj-Pa) 

^ a 

(9) 


where (w) and G°'l'^\uj) are advanced (retarded) 
Green function of boson and fermion, respectively, and 
the subscript iX refers to the enforcement of the local con¬ 
straint on the impurity via Lagrange multiplier A when 
evaluating these correlation functions. The (cu) and 

S j (w) are the self-energies of slave-boson and pseud¬ 
ofermion greater (lesser) Green functions, respectively. 
The NCA expressions for the lesser self-energy of the 
pseudofermion, G^^{uj) = S<(u;)|G^(a;)p, and slave- 
boson, D^{ui) = n<(a;)|iA’’(w)p, are|k [^ [^: 

= J dera{u}-e-fj,a)f{uj-e-Ha)D<{e), 


n''(w) = “X! / Pa)fif^-i^- Pa)G}a{^)- 

^ a d 

( 10 ) 

The physical impurity spectral function, Pcr(<j-’, V), is ob¬ 
tained via the convolution of pseudo-fermion and slave- 
boson Greens function based on Eqs. (I71),®,(inD, and m 
as [13: 

V) = J de[/mT)’'(e)G<(w + e)- 

D<{e)ImG^LJ + e)]. (11) 

The normalization factor Z = -^ J duj[M x (w) — N x 
G^(w)] enforces the constraint, < Q >= 1 with M = 
= 2 here. The current going through the impurity 
therefore reads [13: 


n<(a;) = ;^^ J deTa{e-uj - Ha)f{u} - € +Ha)Gfa{e)- 


liV.T) = I 




ri(w) + r/j(a;) 

X [f{uj + eV/2)- f{u:-eV/2)]. (12) 


( 8 ) 




























FIG. 5: The y/T scaling of equilibrium (linear) conductance 
(G(0,0) — G{0,T))/Br (in unit of 2e/h) versus T/T^ck for 
various r with F = 0.28II, = —0.2D with being 

the two-channel Kondo temperature for r = 0. Here, Br are 
non-universal functions of r. The dashed line is a fit to \/T 
behavior. 


where ra(a;) = r(a; — fia) with a = L,R. The equilib¬ 
rium (linear) conductance is directly obtained via 


G(0,T) 


X 


2 e2 

f 2Td^)T,iu) , 

d/H) 

h 

J TL{uj) + rR{uj)\ 

duj ) 


V = 0). 



(13) 


And the nonlinear conductance G{V) is given by 
Eqs.(6)-(10) form a self-consistent set of Dyson’s 
equations within NCA. We solve these equations self- 
consistently and evaluate Eas. (fTT]) - (ll3D based on these 
solutions. 


III. RESULTS. 

A. Quantum critical point at 2CK-LM phase 
transition: impurity spectral function 

The existence of a quantum critical point (QCP) sepa¬ 
rating 2CK for r < Tc from the LM for r > Vc phases 
exists in the PSG Anderson model has been studied 
extensively [ 2 ^ [ 3 ^. The generic phase diagram is shown 
as in Fig. [TJ We focus here on the transport properties 
for the two-lead setup near criticality both in and out 
of equilibrium, especially on the distinct non-equilibrium 
quantum critical properties (see Fig. [1}. In equilibrium 
(^a = 0), the 2CK-LM phase transition is studied here 
by tuning r of the pseudogap power-law DOS of the leads 
with fixed hybridization parameter Pa and the impurity 
energy Cd- The value of is extracted from the local 
impurity spectral function pa-(oj,V = 0) via solving the 
self-consistent Dyson’s equations. Since €d ^ —U12, the 



FIG. 6: Universal scaling in linear conductance G{0,T/T*) 
(in unit of 2e//i and normalized to T,r{T/T*)'^'^) as a function 
of temperature T/T* near 2GL-LM quantum phase transition 
for various values of r. Parameters are the same as in Fig. [4] 
Here, T*, ctt, and Er are defined in the text. The magenta 
dashed line is a fit to Top Inset: The power-law expo¬ 

nent gt in linear conductance Gqcp{T) close to criticality as 
a function of |r — rc|. Solid line is a fit to a linear relation: 
gt = Pt — arlr — rd with Pt ~ —0.145, ar ~ 1.25. Bot¬ 
tom Inset: Grossover temperature T* versus |r — rc|. The red 
dashed line is a fit to |r — rd^. 


2CK pseudogap Anderson model considered here shows 
particle-hole (ph) asymmetry, giving rise to an overall 
asymmetric impurity spectral function with respect to 
the Fermi energy (see Fig. [2D. For r < Tc, Pa{<^) shows a 
non-Lorentzian Kondo peak, a characteristic of the non- 
Fermi liquid 2CK statejj, [sj. In fact, Pcr(a;) exhibits a 
power-law singularity near w = 0 with an exponent being 
r: Pct(w) ^ a;“’'[2^. With increasing r, the Kondo peak 
gets narrower with reduced spectral weight. For r < Tc, 
however, the Kondo peak splits into two, and the ground 
state is in the LM phase. The critical value Tc ~ 0.115 
for Pa ^ 0.3D, Cd ^ —0.2D (see Fig. |2D. The spectral 
weight of the Kondo peaks gets further suppressed with 
increasing r until it completely disappears. At a finite 
bias voltage for r < Tc, the impurity local DOS shows 
splited Kondo peaks at w = ±U/2 (see Fig.jSD. Note that 
the non-zero local DOS of Pai^J = 0) for r > Tc is due to 
the limitation of the lowest temperature Tq ^ 5 x 10“^D 
we can reach numerically. As T —>■ 0 and for r > Tc, we 
expect a complete dip developed in local DOS such that 

Pa{u! = 0 ) ~ 0 . 
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FIG. 7: The non-equilibrium (non-linear) conductance 
G{y,To) (in unit of 2e//i) versus bias voltage V (in unit of 
D = 1) for various r with F = 0.3D, ed = —0.3D. 


FIG. 8: The 2CK scaling of the non-equilibrium (non-linear) 
conductance (G(0,To) — G{V,To))/Tq'^ (in unit of 2e/ft) ver¬ 
sus (a) [eVIksT)^ and (b) {eV/kBT)°'^ for various r with 
F = 0.3D, td = -0.3D. 


B. Universal scaling in linear (equilibrium) 
conductance near criticality 

More clear signatures of the 2CK-LM quantum phase 
transition can be obtained via linear and non-linear con¬ 
ductances. First, we analyze the linear (equilibrium) 
conductance at finite temperatures but zero bias volt¬ 
age G{V = 0,T). Fig. 0] shows G(0,T) for different r 
with F = 0.28D, /tq = 0, = —0.2Z1, and D = 1. 

For r = 0 it is well known that G(0,T) follows the 2CK 
scaling function^: G(0, T) — G(0, 0) = Br^/TJT^^ for 
T < T 2 CX with ^ 3 X 10“®Z? being the 2CK Kondo 
temperature at r = 0 and being a non-universal con¬ 
stant. Here, we set G(0,0) ~ G(0,ro = 5 x 10“^Z1). 
For 0 < r < Tc, however, the 2CK \/T scaling in 
G(0, 0) — G(0,T) ceases to exist. As shown in Fig. 01 
the ^/T behavior in G(0,0) — G(0,r) is clear for r = 0, 
but it deviates from \/T more with increasing r < Tc- 
We propose that this deviation from the conventional 
2CK behavior for T < T 2 CK for 0 < r < Tc can be due 
to the following two scenarios: (i). the emergence of dis¬ 
tinct universal power-law scaling behaviors when the sys¬ 
tem approaches the 2CK-LM quantum critical regime for 
T > T* with T* being the crossover energy scale above 
which quantum critical behaviors are observed, and (ii). 
the existence of a distinct 2CK scaling form other than 
'/T in conductance in low temperature regime T < T 2 CK 
for 0 < r < Tc with T 2 CK being the two-channel Kondo 
temperature for 0 < r < rc. As indicated in Fig. [5l the 
generic 2CK behaviors for 0 < r < Tc can not be mani¬ 
fested in G(0, 0) — G(0,T) as an universal power-law in 
T as it does for r = 0 case. Therefore, instead of analyz¬ 
ing G(0,0) — G(0, T), we therefore address the above two 
scenarios below via trying the possible scaling behaviors 
of G{T) itself. 


First, as r — > Tc, the existence of a quantum critical 
regime for T > T* requires the divergence of the corre¬ 
lation length ^ in a power-law fashion: ^oc|r — rd”''—>■ 
00 fli- As a result, all thermaldynamical observables, 
including conductances, are expected to exhibit univer¬ 
sal scaling properties. As shown in Fig. 01 we indeed 
find numerically that the linear conductance Gqcp( 0, T) 
shows an universal power-law in T near criticality within 
a temperature range of (approximately) 5 x 10“^Z1 < 
T < 5 X IQ-^D as: 


Gqcp( 0,T) OC = T^T-aTk-r-el^ (-^4) 

where the exponent crp exhibits a linear relation to |r—rd 
with Pt ~ —0.145 and ap ~ 1-25 being non-universal 
constant pre-factors dependent on F and e^, and = 
as r = Tc ~ 0.115. Based on the above analysis, we 
divide G(0,r) by the power-law function T'^t ^ gives: 

G(o,r) = G(0.r)/r- = jjkedd-, 

GQCp(0,r) = GQCp(0,r)/T‘^- (15) 


where Gqcp{0,T) in the quantum critical region be¬ 
comes a constant: Gqcp(0,T) ~ The univer¬ 

sal scaling function G(0, ^) in linear conductance is ob¬ 
tained by rescaling G{0,T) by a non-universal factor Sr 
and T by the crossover energy scale T* (see Fig. El): 


G(0,|^)^ 


G(0, 

(T/r*)'®'r““'rl’—’’cl ’ 


(16) 


where T* is proportional to the inverse of correlation 
length 1/^T OC |r — rd^^, vanishing in a power-law of 
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FIG. 9: Universal scaling in non-linear conductance 

G{V,To)/{V yB r) (in unit of 2ejh) as a function of tempera¬ 
ture VjV* near 2CL-LM quantum phase transition for various 
values of r. Parameters are the same as in Fig. [4] Here, V*, 
CTT, and Br are defined in the text. Black dashed line is a 
power-law fit to the conductance in the 2CK regime. Top 
Inset: The power-law exponent av in non-linear conductance 
Gqcp {V, To) close to criticality as a function of |r — Ccj. Solid 
line is a fit to a linear relation: ay = Pv — av\r — rd with 
Pv ~ 0, ay ~ 2.5. Bottom Inset: Crossover temperature V* 
versus |r — rc|. The red dashed line is a fit to |r — rc|° ®. 


|r —rd with the exponent vt being the correlation length 
exponent: 

T* ox ^ = D\r - (17) 

C.T 


B. Universal scaling in non-linear (non-equilibrium) 
conductance near criticality 

We now add the bias voltage in the leads and study 
the scalings in non-equilibrium conductance near 2CK- 
LM quantum critical point. It is generally expected that 
the scaling behaviors of non-linear conductance near crit¬ 
icality are distinct from that in equilibrium^. The non¬ 
linear conductance G'(y,ro) is obtained at a fixed low 
temperature Tq = 5 x 10“^, symmetrical hybrizidation 
Tl = Tfl = 0.3ZI, and Cd = —0.3D. As shown in Fig. [3 
the non-linear conductance G(V,To) for each value of r 
saturates as U —>■ 0, while a small “kink”-like minima 
in conductance is developed at eV ~ fc^T. For r = 0 
and T < follows the well-known 2CK scaling Q: 

G(0,To) — G(U,To) = where the univer¬ 
sal function behaves as for eV < ksTg, 

and (^ for eV » fesTo. However, for 0 < r < rc as 
the system gets closer to criticality, while behav¬ 

ior remains in G(0,To) — G(U, Tq) for eV ^ ksTo, more 
deviations from the 2CK ( u^To ) ^ behavior are observed 
for ksTfo < eV < kBT 2 CK- Furthermore, similar to the 
case in our analysis for equilibrium conductance, as indi¬ 
cated in Fig. [H the universal 2CK behavior in non-linear 
conductance for 0 < r < rc in general are not mani¬ 
fested in G(0,To) — G(U, To) in the power-law fashion 
as it is for r = 0 case. We therefore analyze the univer¬ 
sal scaling properties based on G(U, Tq) itself rather than 
G(0, To) — G(U, To). We perform the similar analysis here 
to that in equilibrium. As r — >■ r^, we find G(y, To) shows 
a power-law dependence on V in the quantum critical 
regime approximately 10' ^ < V/D < IQ dehned as 
Ggcpiy^To) with power-law exponent linear in |r — rd- 

Gqcp{V,To) OC U'"'" = 


with pt being the correlation length exponent corre¬ 
sponding to the power-law exponent of crossover scale 
T* versus \r — rd- We hnd ~ 4 here, (see Inset of 
Fig. in]). As r —>■ Tc from below, we find a perfect data 
collapse in G(0, ^) ~ const, for 1 < T/T* < 10^, indi¬ 
cating the range of quantum critical region. Surprisingly, 
at low temperatures T < T 2 CK ~ O.IT* where the sys¬ 
tem is governed by the 2CK phase, we find the above 
function G(0, ■^) exhibits as a distinct universal func¬ 
tion power-law scaling function (see Fig. [6|); 

GiO,7^)2CK OC (18) 


with dv ^ otv being non-universal prefactors dependent on 
F and e^. At criticality r = Tc, cry ^ 0, yielding a con¬ 
stant non-linear conductance: Gqcp(H, To)|r,, ~ const.. 
We further perform a re-scaling on G(U, Tq) as: 


G(H,To) 


G(U,To) ^ G(U,To) 


( 20 ) 


where G(V, Tq) becomes a constant in the critical regime. 
We may further dehne the universal scaling function 
G(-^,To) for non-linear conductance via the following 


re-scalings: V V/V*, G(V,To) G(V,To)/Sr (see 
Fig.[9|): 


with <J 2 ck ~ for the parameters we set (or equiv¬ 
alently G{ rp^^ ) OC )<^2ck+<^t'j^ yyg thercfore 

regard this low temperature behavior in G(0, T) as the 
distinct 2CK scaling in equilibrium (linear) conductance 
for the pseudogap 2-channel Anderson model. 


U G(^,To)/S, 

U* ’ ^ (.^)^v-av k-rd 


( 21 ) 


where V* is the crossover energy scale, and a non- 
universal constant pre-factor. Here, we find the exponent 
cry depends linearly on |r — rd with dv ~ 0, ay « 2.5 
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by the best fit of the data (see Inset of Fig. |9l), and 
the crossover scale V* is inversely proportional to the 
correlation length with a power-law dependence on 
the distance to criticality: V* oc oc |r — with 

i^y ~ 0.5 (see Inset of FiglS]). Note that these critical 
exponents out of equilibrium are distinct from those in 
equilibrium, and can be considered as characteristics of 
non-equilibrium quantum criticality. The distinction be¬ 
tween equilibrium and nonequilibrium quantum critical 
properties is expected due to the different role played by 
the temperature and voltage bias near criticality, leading 
to different behaviors in decoherence rate F® (the broad¬ 
ening of impurity DOS) in equilibrium F|n versus out of 
equilibrium Fy [3, [s^ . 

For Tq < V < V*, the system approaches 2CK state 
at a characteristic energy scale V ~ V 2 CK ~ O.IF* be¬ 
low which an universal power-law scaling is observed (see 
Fig. 0: 

Gi^.TohcK cx (22) 

where the exponent ^ 0.23. In the low bias limit 

V < Tq, conductance reaches a constant equilibrium 
value G(0, Tq) and therefore deviates from the 2CK scal¬ 
ing (see Fig. 0. 

Note that we find Eq. for non-equilibrium conduc¬ 
tance to exhibit the same form as that in equilibrium in 
Eq. (IT^ with the same exponent {(J 2 CK = ^ 2 CK ~ 0.23). 
Eurther studies are required to clarify if this relation 
holds in general. Though the exact value (0.23) of the 
above exponents depend in general on physical param¬ 
eters, such as: and F, Eqs. and (|2^ correctly 

reproduce the well-known \/V and \/T behaviors for 
r = 0 in the 2CK regime in G(0,To) — G(E, To) and 
G(0,To) — G(0,T), respectively. 

V. CONCLUSIONS. 

In conclusions, we have studied quantum phase tran¬ 
sitions in and out of equilibrium between two-channel 
Kondo and local moment phases in the two-channel pseu¬ 
dogap Anderson model where the conduction leads show 
a power-law vanishing density of states with exponent 
0 < r < 1. Via Non-Crossing Approximation (NCA), we 
solved self-consistently for the impurity Green’s function, 
and therefore determined the linear and non-linear con¬ 
ductances. The 2CK-LM quantum criticality is reached 
by varying the power-law exponent r of the pseudogap 
conduction bath with fixed lead-dot hybridizations and 
chemical potentials. The linear GiV = 0,r) and non¬ 
linear G(V, To) (To ^ 5 X 10”^T)) conductances show 
distinct universal power-law scalings near criticality with 
different critical exponents. Furthermore, in the 2CK 
regime T,V < T 2 ck^ we also find different characteris¬ 
tic power-law scalings in G(T) and G(y,To) compared 


to the well-known y/T, y/V scalings for G(0,0) — G(0, T) 
and G(V, 0) — G{V, Tq), respectively. Our results provide 
further insights on two-channel Kondo physics and on 
the non-equilibrium quantum criticality in nano-systems. 
Further analytic and numerical investigations are needed 
to address issues on the mechanism behind the different 
scalings between equilibrium and nonequilibrium conduc¬ 
tances. 
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